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Recently, the new strongly intensive measures of fluctuations A and E 
have been proposed. In this publication their properties are tested using 
an example of event-by-event transverse momentum fluctuations. The ob- 
tained values are compared to the long used <i> measure of px fluctuations. 
Several tests are preformed within data produced by fast generators, as well 
as by the UrQMD model. The UrQMD calculations are presented for the 
systems and energies which are planned to be studied in the NA61 /SHINE 
experiment at the CERN Super Proton Synchrotron. 

1. Introduction 

The main motivation of colliding relativistic heavy ions is to create and 
study the properties of the system composed by deconfined quarks and 
gluons (QGP). The data suggest that the energy threshold for deconfine- 
ment (onset of deconfinement) is located at low SPS energies [H [2]. The 
phase diagram of strongly interacting matter can be presented in terms of 
temperature T and baryochemical potential ^b- The bulk of theoretical 
calculations suggest that the phase boundary between hadrons and QGP 
is of first order at large values oi hb, ending in a critical point of second 
order and then turning into a continuous rapid transition at low fi b ■ Lattice 
QCD calculations indicate that the critical point (CP) can be located in the 
SPS energy range, i.e. T"^^ = 162 ± 2 MeV, /ig^ = 360 ± 40 MeV [3] or 
{T^^,^l%^) = (0.927(5)Tc,2.60(8)rc) = (~ 157, ~ 441) MeV gj, where 
is the critical temperature of hadron gas -f-)- QGP transition at vanishing 
baryochemical potential. 

The analysis of dynamical fluctuations can be an important tool for lo- 
calizing the phase boundary and the critical point. In particular, significant 
transverse momentum and multiplicity fluctuations are expected to appear 
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for systems freezing-out close to CP [5j. The position of the freeze-out point 
in the phase diagram can be moved by varying the cohision energy and the 
size of the cohiding nuclei |6]. A non-monotonic evolution of fluctuations 
with such parameters can serve as a signature for the phase transition and 
the critical point. 

In fact, these considerations motivated an extensive program of fluctu- 
ation studies at the SPS and RHIC accelerators. The NA49 [7] experiment 
reported non-monotonic behavior of average px and multiplicity fluctuations 
at the top SPS energy [8j. This intriguing result might be a first hint of the 
critical point. Therefore the efforts to look for CP will be continued within 
the NA61/SHINE [9] project where a 2D (energy and system size) scan 
of the phase diagram will be performed. Figure [T] shows the hypothetical 
chemical freeze-out points in the NA61 experiment. 

Dynamical fluctuations can 
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Fig. 1. Positions of chemical freeze-out points 
obtained within hadron gas model by fitting 
NA49 data |5] (blue squares). Circles are those 
expected in NA61. Taken from [lOj . 



be measured by use of event- 
by-event methods. However, 
while measuring event-by-event 
fluctuations in nucleus-|-nucleus 
(^-|-A) collisions one should re- 
member about a trivial source 
of fluctuations caused by event- 
by-event changes of the col- 
lision geometry. Therefore, 
a suitable choice of statistical 
tools for the study of event-by- 
event fluctuations is really im- 
portant. In Ref. [TI] a strongly 
intensive measure $ was in- 
troduced. In a superposition 

model ^ does not depend on the number of "sources" (e.g. wounded nucle- 
ons in the Wounded Nucleon Model [12] ) composing ^ + j4 collision (inten- 
sive measure) and on the fluctuations of this number of sources (strongly 
intensive measure). In addition, in thermodynamical models ^ does not 
depend on volume and volume fluctuations provided that temperature and 
chemical potential are constant. The ^ measure was already used by the 
NA49 experiment to calculate transverse momentum fluctuations {^p^) 
[131 E] • Recently, two new classes of strongly intensive measures have been 
proposed: A and S [15j. In fact, previously proposed $ belongs to S-family 
measures. In heavy ion experiments, the use of strongly intensive measures 
of fluctuations, such as A and S, can be a remedy for an imperfect cen- 
trality selection of ^ + ^ collisions. Therefore, the NA49 experiment was 
able to use relatively wide centrality bins while studying measure (up 
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to 0-15% centrality) , whereas the analysis of multiphcity fluctuations had 
to be hmited to 1% most central interactions only [8]. 

In this paper first basic tests of the newly proposed A and T, measures 
win be presented for transverse momentum fluctuations. The obtained val- 
ues will be compared to the long used ^ measure oipx fluctuations. Several 
effects have been studied for events generated by use of so-called fast gen- 
erators. Moreover, the analysis within a much more complex UrQMD3.3 
model |16) will be shown. The UrQMD calculations have been done for 
the systems and energies which are planned to be studied in the CERN 
NA61/SHINE experiment (see Fig. d]). 

2. Strongly intensive measures 

We call intensive quantities those ones which do not depend on the 
volume of the system. In contrary, extensive quantities (for example mean 
multiplicity or variance of multiplicity distribution) are proportional to the 
system volume. Note, that it is useful to extend the notion of intensive and 
extensive quantities to the Wounded Nucleon Model. Namely, the intensive 
quantities can be called those ones which are independent of the number of 
wounded nucleons, and extensive ones those which are proportional to the 
number of wounded nucleons. The ratio of two extensive quantities is an 
intensive quantity [15J. Therefore the ratio of mean multiplicities, as well 
as the scaled variance of multiplicity distribution uj = {{N^) - {N)^)/{N), 
are intensive measures. In fact, due to its intensity property lv measure 
is quite commonly used to determine multiplicity fluctuations in heavy ion 
experiments. 

There is one more important problem which one should not forget about. 
In high energy heavy ion collisions the volume of the produced matter cannot 
be fixed. In fact the system volume changes significantly from event to event. 
Therefore it is very important to be able to measure the properties of the 
created matter independently of its volume fiuctuations. The quantities 
which allow this are called strongly intensive measures. They do not depend 
on the volume and on volume fluctuations. 

Mean multiplicities ratios are both intensive and strongly intensive mea- 
sures. The situation is however more difficult for fluctuation analysis. The 
scaled variance of multiplicity distribution is an intensive measure but not 
strongly intensive. Quite long ago a strongly intensive measure <I> was first 
introduced [11]. In the recent paper [15] is was shown that there are at 
least two families of strongly intensive measures: A and S. The previously 
known $ measure belongs to S-type family. They can be calculated for any 
two extensive quantities. In this paper A, S, and ^> measures calculated for 
particle multiplicity, N, and sum of their transverse momenta modules will 
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be tested. 



2.1. <5pj, measure 

The $ measure m was already successfully used by NA49 to determine 
transverse momentum fluctuations ('^py) [131 I14j . Following the authors 
of [11] one defines the single-particle variable Zpj, = pr —pf with the bar 
denoting averaging over the single-particle inclusive distribution. As seen 

= 0. Further, one introduces the event variable , which is a multi- 
particle analog of Zp^, defined as Zp^ = J2f=iiPTi —pf), where the summa- 
tion runs over particles in a given event. Note, that (Zp^) = 0, where (...) 
represents averaging over events. Finally, the $pj, measure is defined as 



is a strongly intensive measure and therefore \i A + A collision is 
represented by an incoherent superposition of independent nucleon-l-nucleon 
(A^-|- A^) interactions (superposition model), then $py has a constant value, 
the same for A + A and N + N interactions. This implies that, in particular, 
$pj, does not depend on the impact parameter (centrality) , if the A + A 
collision is a simple superposition oi N + N interactions. 

Another property of this measure is that $pj, vanishes when the sys- 
tem consists of particles that are emitted independently (no inter-particle 
correlations) and the single particle momentum spectrum is independent of 
multiplicity. In contrary, A and S measures, which will be shown below, do 
not assume zero values for independent particle production. 

2.2. and Ti^^ measures 

Let A and B be two fluctuating extensive quantities. Then A^^ and 
S"^^ can be defined 1151: 



A^"" = {C)-'[{B)iOA-{A)iOB], (2) 

E^^ = {C)-'[{B)u;a + {A)u;b - 2{{AB) - {A){B))], (3) 

where: 

= jj^ (4) 

and 

= (5) 
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are the scaled variances of two fluctuating extensive quantities A and B. 
(C) is the average of any extensive quantity e.g., {A) or (B). 

There is an important difference between A^^ and T,^^ . Only the first 
two moments: (A), {B), and (^^), (-B^) are required to calculate A"^^, 
whereas includes the correlation term (AB) — {A){B). Thus and 
^A_B jxieasures can be sensitive to several physics effects in different ways. In 
publication [TS] all strongly intensive quantities including correlation term 
are named the S family, and those including only mean values and vari- 
ances the A family. The already used $ measure belongs to S-type family. 
The definitions of new quantities A^^ and Y^"^^ are, however, more general 
because one can use here any two extensive measures A and B, whereas in 
the original definition of ^ one of them was fixed to be particle multiplicity. 
Nevertheless, in this paper one of the extensive measures will be again mul- 
tiplicity, in order to compare the results with known measure <I>. Finally, 
one should mention that A'^^ and T,^^ have also different properties with 
respect to exchange A and B: S^-^ = S^^ and A^^^ = -A^^^ 

For the analysis of transverse momentum fluctuations one can use |15j : 
A = X = J2iLi (where Xi = pT,i and the summation runs over all accepted 
particles in a given event), B = N, C = N, and X2 = J2iLi ^i- Then we 
obtain: 



(iV) 

and 



(N) 




(6) 



= J^m^x + {X)UN - 2{{XN) - {X){N))]. (7) 
The Tj^^ measure can be also expressed [15j using known <I>pj, quantity: 

W. (8) 



(X) {X) (N) 



In this paper S is calculated using Eqs. ([8]). The statistical errors on 
$pj,, A"'^^ and T^-^^ were estimated as follows. The whole sample of events 



was divided into 30 independent sub-samples. The values of '^py, A , and 
were evaluated for each sub-sample and the dispersions {D^, Da, and 
-Ce) of the results were then calculated. The statistical error of ^pj, (A^^ 
or S^^) is taken to be equal to D^/^/m {Da/V^ or 
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3. Results of fast generators 

3.1. Independent particle production 

The basic properties of <^p^ , A^-^ and T,^^ were tested by use of Monte 
Carlo models (so-called fast generators). Each interaction (event) was com- 
posed by a given number of identical single sources. For each source the 
number of particles was generated from the Poisson distribution with a 
mean value of 5. The particle transverse momentum was generated from 
exponential transverse mass spectrum with inverse slope parameter ("tem- 
perature") T = 150 MeV. The number of sources composing an event was 
either constant (circles in Fig. [2]) or selected from Poisson (triangles) or 
from Negative Binomial distribution (squares). For Negative Binomial dis- 
tribution its dispersion \JVar{Ns) was large and taken to be equal {Ns)/2. 

Figure[2]presents ^pj,, A"^^ and Tj-^^ measures versus the number or the 
mean number of sources composing one event. As the simulation performed 
here represent independent particle production measure is consistent 
with zero. In contrary, A'''-'^ and have non-zero values because their 

definitions do not assume zero for independent particle production. The 
circles nicely confirm the intensity property of all three measures, whereas 
the triangles and the squares show that these quantities are also strongly 
intensive (do not depend on the system size and the system size fluctua- 
tions). One can also mention here that for a constant number of sources 
per event (circles in Fig. [2]) the scaled variance of multiplicity distribution 
ui = {{N'^) — (N)'^) / (N) is close to 1 in the whole range of the horizontal 
axis. For Poisson number of sources distribution ui is approximately 6 also 
in the whole range of the mean number of sources per event. For Negative 
Binomial distribution of the number of sources ui increases from about 7 (for 
on average 5 sources per event), through 126 (for on average 100 sources 
per event) up to approximately 1000 (for on average 800 sources per event). 
Figure [2] shows that within statistical errors (which can be sometimes large) 
all three measures are strongly intensive even if the multiplicity distribution 
is extremely and abnormally wide. 

3.2. "Temperature" fluctuations 

In the next simulation for each single source the number of particles was 
again selected from the Poisson distribution with a mean value of 5. The 
particle transverse momentum was generated from exponential transverse 
mass spectrum with average inverse slope parameter (T) = 150 MeV. The 
T parameter was generated separately for each single source (source-by- 
source T fluctuations) from Gaussian shape with dispersion cjt = 25 MeV. 
Finally, the number of sources composing an event was generated from the 
Poisson distribution. The results are presented in Fig. O The effect of 
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Fig. 2. $p^, A^^^ and T,^'" versus number or mean number of sources composing 
one event. Example of independent particle production. 



"temperature" fluctuations results in positive ^pj, values and in higher (then 
in Fig. [2|) A''^^ and values. As the sources are identical (superposition 
model) ^pj,, A^'^ and Tj-^^ measures do not depend on the mean number 
of sources composing an event. Therefore Fig. [3] also confirms that all three 
measures are strongly intensive. 
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Fig. 3. $pri A and S versus mean number of sources composing one event. 
Simulation of the effect of source-by-source inverse slope parameter (T) fluctua- 
tions. 



In the next step source-by-source T fluctuations from the previous sim- 
ulation were replaced by event-by-event T fluctuations (T parameter was 
the same for all sources composing a given event but was varied between 
events). Here the T parameter was generated separately for each event 
from Gaussian shape with fj-r = 25 MeV. Again, the number of sources 
composing an event was generated from the Poisson distribution. The re- 
sults are shown in Fig. HI In this simulation the values of ^pj,, A-^'^ and 
j^XN i22crease with increasing the number of sources composing an event. 
It suggests that in real experimental data the effect of event-by-event " tem- 
perature" fluctuations should be better detectable for more central (or for 
heavier A) collisions. One should also mention here that the relationship 
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between temperature and multiplicity (or volume) fluctuations was studied 

in [m m HI] . 
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Fig. 4. A and S versus mean number of sources composing one event. 

Simulation of the effect of event-by-event inverse slope parameter (T) fluctuations. 

In the last simulation different widths of T parameter distribution were 
considered. As usually, for each single source the number of particles was 
selected from Poisson with a mean value of 5. The particle transverse mo- 
mentum was generated from exponential transverse mass shape with average 
inverse slope parameter (T) = 150 MeV. The T parameter was varied from 
event to event following Gaussian distribution with dispersion ut- In order 
to avoid negative T values only events within T = 150 it Sax MeV were 
accepted. Finally, the number of sources composing an event was gener- 
ated from the Poisson distribution with a mean value of 100. The results 
are shown in Fig. [5j As expected, the values of all three fluctuation mea- 
sures increase when event-by-event "temperature" fluctuations are stronger 
(higher ax)- 
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Fig. 5. and E'^^ versus width of T parameter distribution. 
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The influence of "temperature" fluctuations 
on transverse momentum fluctuations was al- 
ready studied in [20]. Following the authors of 
[20j one can easily derive that event-by-event 
"temperature" fluctuations result in <I>pj, equal 

64 + 2(r)2. 



4(2 + 4(Ar2)/(Ar)) + 2(r)2 - 
For the scaled variance of multiplicity distribu- 
tion u) ^ 6 (taken from Monte Carlo simulations 
in Fig. [5]) the formula can be rewritten as: ^pj, = 
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Fig. 6. versus relative 
width of T parameter dis- 

24(13 + 2{N)) + 2(T)2 - ^^64 + 2(T)2. This tribution. See text for de- 
function, for (T) = 150 MeV and (N) = 500, is tails, 
drawn solid line in Fig. [6l together with re- 
sults from Monte Carlo simulations (taken from 

Fig. [5] left). The Monte Carlo results are in a very good agreement 
with the analytical formula. In publication [20] it is also suggested that 
for sufficiently small T parameter fluctuations, it is when {N){T)'^ ^ 
(iV2)((r2) - (r)2), and for Poisson (w = 1) multiphcity distribution 

simplifies to the formula y/2{N) ^'^ (y)"^^ ■ This simplified function, for 
(T) = 150 MeV and (A^) = 500, is drawn as a dashed line in Fig. [6l 
As seen, the simplified function is indeed valid only for small T parameter 
fluctuations and starts to deviate from Monte Carlo simulations for ax/ {T) 
higher than approximately 3 percent. 



4. Results of the UrQMD model 

In order to calculate A^^, S^^ and ^p^ the UrQMD3.3 model [H] 
was used. The model was already used to study event-by-event fluctua- 
tions of average transverse momentum |2H [22] , charged particle ratio |23] , 
and particle number [24J.The UrQMD generator is a microscopic transport 
model producing hadrons via formation, decay, and rescattering of reso- 
nances and strings. The UrQMD approach simulates multiple interactions 
of both target/beam nucleons and newly produced particles, excitation, and 
fragmentation of color strings and the formation and decay of hadronic reso- 
nances. The samples of minimum bias events were generated for the systems 
which will be analyzed in the NA61/SHINE experiment, namely Be + Be, 
Ar + Ca, Xe + La at beam energies 13^1, 20A, 30A, 40^, 80A, and 158^ 
GeV (Fig. [1]). These beam energies correspond to center-of-mass energies 
for nucleon-|-nucleon pair y^s/vjv equal 5.12, 6.27, 7.62, 8.76, 12.3, and 17.3 
GeV, respectively. Only inelastic A + A collisions were considered, it is the 
generated events with the number of collisions equal zero and events with 
the number of inelastic collisions equal zero were removed from the sample. 



10 



The remaining events which had the final state multiplicities equal to the 
sum of the nucleons in both colliding nuclei were also rejected. 

Prom such prepared minimum bias events only 20% of the most central 
collisions were selected (in agreement to NA61 data taking plans). This 
centrality selection corresponds to impact parameter cuts in UrQMD: b < 
1.96 fm for Be + Be, b < 3.81 fm for Ar + Ca, and b < 5.71 fm for Xe + La. 



4-1. Influence of high-pT particles (at SPS energies) 

Figure [3 presents A"'^^ and T,-^^ measures obtained for all charged 
particles produced in UrQMD. Open points are those without any kinematic 
restrictions. Full points correspond to results with 0.005 < px < 1-5 GeV/c 
cut. Such cut was applied by the NA49 experiment |T3l[lll[25] in order to 
eliminate the possible effect originating from hard interactions (the lower 
cut was due to momentum resolution of the detector). Figure [8] shows mean 
multiplicities of charged particles with and without 0.005 < pt < 1-5 GeV/c 
cut. Although at SPS energies there is only a small fraction of particles 
coming from hard interactions (compare open and full points in Fig. [8| 
such high-pT^ hadrons can increase the values of all considered fluctuations 
measures (Fig. [7|). For Xe + La collisions at the SPS this increase is on the 
level of 15-19% for all three fluctuation measures. 
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Fig. 7. $prj A and S measured for all charged particles produced in 20% 
most central A + A collisions in UrQMD. 



In Fig. [7] one also observes a slight decrease of ^pj, and T,-^^ with 
increasing energy, while, in contrary, A'''"^ increases. It is also worth to 
notice that the values of all three measures are lower for lighter systems 
(increase when going from Be + Be to Xe + La). In a simple superposition 
model no system size dependence is expected. However, the increase of ^pj,, 
A-^^ and T^-^^ for heavier systems may be due to event-by-event impact 
parameter fluctuations. The sample of 0-20% central collisions may be not 
narrow enough (especially for the heaviest Xe + La) and the characteristics 
of more and less central collisions may be significantly different leading 
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Fig. 8. Mean multiplicities of charged particles produced in 20% most central A + A 
collisions in UrQMD. 



to increased event-by-event fluctuations. In the next sections it will be 
shown that indeed the values of fluctuation measures are smaller when the 
centrality is restricted to 7.2% most central data. 

^.2. Influence of hadrons different than vr, K, anti{p) 

In the event-by-event analysis of experimental data [131 [3 ES] typi- 
cally only charged particles originating from the main vertex are used. It 
practically means that only main vertex pions, protons, kaons and their an- 
tiparticles are used in the analysis, because particles coming from the decays 
of for example A, S, H and are believed to be rejected by a set of track 
selection criteria. The time scale of the simulations performed within the 
UrQMD model did not allow for weak decays, therefore the UrQMD anal- 
ysis of charged pions, protons, kaons and their antiparticles can be directly 
compared to the analysis of experimental data. 

Figure M compares ^p^, , A^^ and T,^^ calculated for all charged parti- 
cles produced in UrQMD (full symbols) and for only charged pions, protons, 
kaons and their antiparticles (open symbols). The small difference between 
these two cases is a reflection of the fact that the majority of final state par- 
ticles produced in relativistic heavy ion collisions are pions, (anti)protons 
and kaons. 



4-3. ^pj, - dependence on energy and charge combination 

Figure [TOl shows the energy dependence of ^p^, for different charge com- 
binations: all charged particles, negatively charged and positively charged. 
The open symbols in the right panel represent positively charged particles 
where protons were removed from the sample. As already discussed in Fig. 
[7] a decrease of ^p^, with energy can be observed for all charged particles. 
This decrease is even more pronounced when looking at positively charged 
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Fig. 9. ^pT, A"^^ and E'^^ measured for all charged particles and for charged 
pions, kaons and (anti)protons produced in 20% most central A + A collisions in 
UrQMD. 



particles (Fig. [TOl right). In a contrary, negatively charged particles show 
values close to zero and only a very small increase with energy (see inset of 
Fig. [TU] middle) can be observed for heavier systems. The most important 
difference between positively charged and negatively charged particles is the 
presence of protons in the sample of positively charged. Therefore the $pj, 
values were calculated also for positively charged particles without protons 
(open points in the right panel). Their values are consistent with those 
for negatively charged particles thus confirming that the effect indeed orig- 
inates from protons. It is also worth to remind that the highest effect from 
protons is observed for the lowest energies where the fraction of protons is 
the highest. For Xe + La UrQMD data the fraction of protons is about 
35% of all charged for 13A GeV and 15% of all charged for 158A GeV. For 
Ar + Ca data these numbers are 36% and 18% for 13^ GeV and 158yl GeV, 
respectively, and for Be + Be collisions 48% and 24%. 
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Fig. 10. Energy dependence of for different charge combinations of particles 
produced in 20% most central A + A collisions in UrQMD. 

In the next set of plots (Fig. Ilip the ^p^, values are shown for forward 
rapidity only (1.1 < y* < 2.6, where y* is the particle rapidity calculated 
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in the center-of-mass reference system assuming pion mass). The same 
rapidity cut was used by the NA49 experiment in the analysis of transverse 
momentum [T^ and multiphcity [25] fluctuations. Additionally, only 
particles with y* < y^g^^ — 0.5 were accepted iyleam rapidity of the 

beam calculated in the center-of-mass reference system). This cut allows to 
get rid of the effect of event-by-event impact parameter fluctuations when 
restricting the analysis to forward rapidity only. The details of this cut 
were discussed in [22\ and in particular it was shown that event-by-event 
fluctuations of the number of protons in the forward hemisphere and the 
number of protons that are closer to the production region can lead to 
increased <I>r,„ values. 
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Fig. 11. Energy dependence of at forward rapidity for different charge com- 
binations of particles produced in 20% most central A + A collisions in UrQMD. 
Additional cut y* < j/beam ~ ^■'^ applied. 
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Fig. 12. Mean multiplicities of all charged, negatively charged and positively 
charged particles produced at forward rapidity in 20% most central A + A col- 
lisions in UrQMD.. Additional cut y* < yleam ~ 0-5 was applied. 

Figure [11] shows that ^p^, values measured at forward rapidity are much 
smaller than those for complete rapidity region (see Fig. [10] for a com- 
parison). Also the mean multiplicities are significantly smaller (Fig. I12|) . 
The negatively charged particles show <I>pj, consistent with zero for all three 
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studied systems (the slight increase with energy observed for complete ra- 
pidity region is not seen here any more). The positively charged particles 
also result in $pj, values close to zero and a weak increase with increasing 
energy can be observed for Xe + La system only. 

4.3.1. Centrality restrictions 

Figure [TO] suggests that for 20% most central interactions the increased 
values at lower SPS energies are mainly due to the significant fraction of 
protons present in the samples. The centrality bin 0-20% is a relatively wide 
one and therefore the observed increase of may be due to event-by-event 
impact parameter fluctuations (corresponding to event-by-event fluctuations 
in the number of participating protons). One can suppress this effect by 
selection of narrower centrality bins. In Fig. [13] the centrality of Xe + La 
at 13^ GeV beam energy is restricted from 0-20% (rightmost points) down 
to 0-1% most central (leftmost points). As the same (0-20%) event sample 
was used points may be correlated. The negatively charged particles do 
not show any dependence on a/atotai- As expected, the values for all 
charged paricles and for positively charged ones decrease when going to 
more central collisions, reaching the values similar to that for neagatively 
charged at aproximatelly 7% most central interactions. 
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Fig. 13. ^pj. and mean multiplicities as fmiction of percent of total inelastic cross 
section for Xe + La collisions at the lowest SPS energy. Note: the values and their 
errors are correlated. 

Figure [13] confirms that the centrality range (0-20%) of the events planned 
to be recorded by the NA61 experiment may need futher restrictions in or- 
der to calculate transverse momentum fluctuations. Therefore in Fig. [TJ] 
the same results as in Fig. [10] are presented but additional cut on selection 
of 7.2% most central collisions was applied. This cut corresponds to cen- 
trality selection done by NA49 in the study of the energy dependence of ^p^ 
for Pb+Pb interactions [H]. In Fig. [H] additional points for UrQMDl.3 



^ In UrQMD simulations centrality 0-7.2% corresponds to impact parametr cut b < 1.17 
fm for Be + Be, b < 2.29 fm for Ar + Ca, and b < 3.42 fm for Xe + La interactions. 



15 



Pb+Pb collisions are included. For 7.2% most central interactions ^pj, mea- 
sure is close to zero but only for negatively charged particles. For positively 
charged it is close to zero but only for heavier systems {Pb + Pb and 
Xe + La), however for light systems (especially Be + Be) we still observe 
a significant increase at lower SPS energies. This increase is not visible 
any more when protons are removed from the sample. As the centrality is 
already restricted here we can guess that there may be another yet source 
of correlations (conservation laws?) that produces positive ^p^, values for 
central interactions at lower energies. This hypotesis needs further investi- 
gations however it is supported but the recent observation [26| that also in 
UrQMD p + p interactions (no centrality restrictions required) ^pj, values 
for all charged and for positively charged particles are increased at lower 
SPS energies, wherease the negatively charged ones are consistent with zero 
in the whole SPS energy range. 
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Fig. 14. Energy dependence of for different charge combinations of particles 
produced in UrQMD. 7.2% most central interactions are selected. 



4-4- A - dependence on energy and charge combination 

In Fig. [15] the energy dependence of A'''-^ is shown for 20% most central 
A + A interactions. Three different charge combinations are included: all 
charged particles, negatively charged and positively charged. The open sym- 
bols in the right panel represent positively charged particles where protons 
were removed from the sample. In contrary to ^pj, measure (see Fig. [TO]l 
the A"'*''^ measure shows monotonic increase with increasing energy for all 
three charge combinations and for all studied systems. Similarly to the 
values of A"'*'^ for positively charged particles are comparable with those 
for negatively charged ones provided that protons are removed from the 
sample. 

In the next set of plots (Fig. [T6|l the A"'^^ values are show for forward 
rapidity only (1.1 < y* < 2.6). Additionally, only particles with y* < 



Very weak increase of $pj, with energy may be observed for Pb + Pb and Xe + La. 



16 



> 3000 - 0.005 < < 1 .5 GeV/c 
" All charged 



• Be+Be 
■ Ar+Ca 
A Xe+La 



^ [GeV] 



- 0.005 <p <1.5GeV/c600 




Negatively charged 




200 


■e : ' 


• Be+Be 




5 10 15 


■ Ar+Ca 


A Xe+La 




• • " : 
r ' ' * ' 


I 


5 10 


15 20 
\|s^ [GeV] 



0.005 <p <1.5GeV/c 


• Be+Be 


Positively ciiarged 




■ Ar+Ca 


□pen - without proto 




A Xe+La 














, . - ' 






. ■ ■ ■ 

• ? ? • 


■ 


8 




□ 






\|s^ [GeV] 



Fig. 15. Energy dependence of A^''^ for different charge combinations of particles 
produced in 20% most central A + A collisions in UrQMD. 



ybeam~^-^ Were accepted. As seen A values measured at forward rapidity 
are much smaller than those for complete rapidity region (see Fig. [15] for 
a comparison). The negatively charged particles show only a very weak 
increase and then saturation of with increasing energy. Positively 

charged particles at forward rapidity show decrease of A'''"^ values with 
energy for heavier systems {Ar + Ca, Xe + La). This is in contrary to 
(see Fig. [TT]) . which showed an increase with energy for positively charged 
particles in Xe + La interactions. 
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Fig. 16. Energy dependence of A"^^ at forward rapidity for different charge com- 
binations of particles produced in 20% most central A + A collisions in UrQMD. 
Additional cut y* < ylf,cmi ~ 0-^ applied. 



Similarly to what was done for ^pj, the centrality was also restricted 
when calculating A"''^^. In Fig. [T7]the centrality of Xe + La at 13^ GeV 
beam energy is restricted from 0-20% (rightmost points) down to 0-1% most 
central (leftmost points). The negatively charged partciles show a small in- 
crease of A'''-^ with increasing a/atotai {^pt values were constant - see Fig. 
[T3|) . The A"**-^ values for all charged paricles and for positively charged ones 
decrease when going to more central collisions, thus confirming that a signif- 
icant contribution to A"'*-^ may be indeed from event-by-event fluctuations 
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of the number of participating protons. 

Figure [18] presents the energy dependence of 
^XN £qj, 7_2% most central ^4+^ colhsions (addi- 
tional points for UrQMDl.3 Pb+Pb collsions are 
included). In the case of Xe+La and Ar+Ca data 
the 7.2% most central interactions result in A"''-^ 
values much smaller than those for 20% most cen- 
tral (see Fig. [15] for a comparison). For Be + Be 
interactions such drop is much less pronounces 
and in fact the magitudes for 7.2% and 20% most 
central interactions are comparable. For all sys- 
tems the values of A"'*"^ for 7.2% most central 
collisions are higher for positivey charged part- 
ciles than for negatively charged ones and the re- 
jection of protons can bring the magintudes of 
positively charged ones close to these for nega- 
tively charged. In contrary to what was observed 
for for very central collisions A^^ for neg- 
atively charged particles show a rudimentary en- 
ergy dependence. All systems and charge combinations show slight increase 
of A'''-''^ with energy and only central Pb-|-Pb simulations seem to exhibit a 
very weak decrease/plateau of A'''"^ with energy. The different behaviour 
of Pb + Pb collisions may be connected with the fact that the earlier ver- 
sion UrQMDl.3 was used for Pb + Pb simulations, wherease UrQMD3.3 was 
employed to simulate Be + Be, Ar + Ca and Xe + La interactions. 



Fig. 17. as function 

of percent of total inelas- 
tic cross section for Xe + 
La collisions at the low- 
est SPS energy. Note: the 
values and their errors are 
correlated. 
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Fig. 18. Energy dependence of A'^^ for different charge combinations of particles 
produced in UrQMD. 7.2% most central interactions are selected. 



In a simple superposition model all three fluctuation measures A 



XN 



^XN 



and ^pj, should have the same value, independent of the system size 
(the number of sources). Therefore when canceling trivial effect of event- 
by-event impact parameter fluctuations (restricting to 7.2% most central 
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interactions) we should expect the same values of fluctuation measures for 
all studied systems. When considering "most save" sample of negatively 
charged particles one can see that it is indeed true for measure (see 
middle panel of Fig. [T^ . whereas for A"''^^ (middle panel of Fig. [T8|) is 
is true only for lower SPS energies. As the results from the fast genera- 
tors in the previous sections showed that A"''^^, and are indeed 
strongly intensive measures, it may suggest that the UrQMD model intro- 
duces somehow a small but systematic deviation from a simple superposition 
model beginning from middle SPS energies. <&pj, measure can be less sensi- 
tive than A""*"^ to such deviation and therefore its values are the same for 
all systems in the whole SPS energy range. 

In the models the deviation from a simple superposition model (Wounded 
Nucleon Model) may be due to changes of the shape of the single particle 
spectrum and/or due to changes of fluctuations. It was suggested [27j that 
an important difference between $ and A or S is that in case of $ measure 
the single particle term is subtracted (the second term of <1> definition in Eq. 
[T]), thus making $ sensitive only to changes of fluctuations and insensitive 
to changes of the shape of the single particle spectrum. Such subtraction 
is not done for A and S, thus in this respect they are equivalent to the 
flrst (fluctuation) term of As a consequence for an uncorrelated particle 
production $ is zero and independent of a shape of the single particle spec- 
trum. In a contrary, the A and S are non-zero and their values depend on 
the shape of the single particle spectrum. If the UrQMD model introduces 
deviations from a simple superposition model by modifying the shape of the 
single particle spectrum only, then the large deviation from a superposition 
behavior should be visible for A and S only but not for <I>. Therefore, in 
the analysis of experimental data a simultaneous measurement of all three 
quantities A, and E may help to understand the origin of superposition 
model violation. 



4.5. S - dependence on energy and charge combination 

Figure [19] shows the energy dependence of for 20% most central 

A + A interactions. Three different charge combinations are presented: all 
charged particles, negatively charged and positively charged. The open 
symbols in the right panel represent positively charged particles without 
protons. The samples of all charged and positively charged partciles show 
a monotonic decrease of S"^^ with increasing energy, whereas negatively 
charged rise and reach saturation at middle SPS energies. Similarly to 
and A"'*-^ the values of for positively charged particles became 

comparable with those for negatively charged ones provided that protons 
are removed from the sample. 
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Fig. 19. Energy dependence of S^''^ for different charge combinations of particles 
produced in 20% most central A + A collisions in UrQMD. 



In Fig. [50] the values of Ti^^ are shown for forward rapidity only (1.1 < 
y* < 2.6). Additionally, only particles with y* < yl^^^ — 0.5 were accepted. 
T;^^ values measured at forward rapidity are slightly smaller than those 
for complete rapidity region. All three charge combinations and all systems 
show an increase and then saturation of with increasing energy. 
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Fig. 20. Energy dependence of T,^^ at forward rapidity for different charge com- 
binations of particles produced in 20% most central A + A collisions in UrQMD. 
Additional cut y* < yl^^am ~ 0.5 was applied. 

Figure ED presents the values of T,-^^ when the centrality of Xe + La 
at 13A GeV beam energy is restricted from 0-20% (rightmost points) down 
to 0-1% most central (leftmost points). The negatively charged particles 
exhibit nearly flat dependence of on a/atotai (similarly to $pj, - see 

Fig. [T3|) . The values of T;^^ for all charged paricles and for positively 
charged ones decrease when going to more central collisions and reach a 
plateau for approximately 10% most central collisions. There is however 
one significant difference between the centrality dependence of or A'''-^ 
and Ti^^ for Xe + La at 13^ GeV. For extremely central collisions the 
values of both and A'''-^ are similar for positively, negatively and all 
charged particles (see Figs. [T3] and [T7|) . In contrary, the values of T,-^^ 
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for very central collisions are much different for negatively charged particles 
than for positively charged and all charged. It might indicate the presence of 
additional source of correlations in the sample of positively charged particles, 
which was not detected by use of ^p^, and A"'''^ measures. The origin of 
such correlation(s) was not investigated here. 

The last set of plots (Fig. [22]) presents the 
energy dependence of for 7.2% most cen- 

tral A + A collisions. For all systems the values 
of Tj-^^ are higher for positivey charged partciles 
than for negatively charged ones and, as usually, 
the rejection of protons brings the magintudes of 
positively charged ones close to these for nega- 
tively charged. 

The 7.2% most central collisions showed a dif- 
ference between the values of fluctuation measures 
for positively charged particles and for positively 
charged without protons. In case of <&pj, (Fig. [HI 
this difference is significant only for light systems 
{Be + Be, Ar + Ca) and in case of A^^ and 
such deviation is observed for all systems. There- 
fore, the results suggest that either the centrality 
bin 0-7.2% is still not narrow enough (it may be 
the case for A'''"'^, see Fig. [T7|l and event-by-event impact parameter (num- 
ber of protons) fluctuations can be still present or there is another yet source 
of correlations visible for positively charge particles only (then protons, of 
course, would have a significant contribution to such correlation). The T,^'^ 
measure can be much more sensitive to such kind of correlation (conservation 
laws?) because a difference between T,-^^ for positively charged particles 
and for positively charged without protons is visible for all systems. In case 
of ^prp, as already mentioned, only smaller systems show this deviation. 
The origin of such possible correlation is still under study and in particular 
the $pj, and T,^'^ measures will be carefully analysed for the energy scan 
of UrQMD p + p cohisions [2^] . 

Finally, in Fig. [22} one observes that for 7.2% most central collisions 
the T,-^'^ values for negatively charged particles are nearly independent of 
energy. There is however rather significant system size (A) dependence 
suggesting the deviation of UrQMD from a simple superposition model. 
Such deviation was also observed for A^^ measure (middle panel of Fig. 
\TE\i . however only at higher SPS energies. 
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Fig. 21. Y.^^ as function 
of percent of total inelas- 
tic cross section for Xe + 
La collisions at the low- 
est SPS energy. Note: the 
values and their errors are 
correlated. 
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Fig. 22. Energy dependence of S'''^^ for different charge combinations of particles 
produced in UrQMD. 7.2% most central interactions are selected. 



5. Summary and conclusions 

In this publication the recently proposed A"''-^ and Tj-^^ measures of 
transverse momentum fluctuations are tested and compared to ^p^, mea- 
sure. The fast generator simulations confirm that all three measures are 
both intensive and strongly intensive. The results suggest that all three 
measures are also very sensitive to the effect of event-by-event "tempera- 
ture" fluctuations. Such effect should be better visible for heavier systems. 

Some preliminary tests were performed also within the UrQMD model, 
which showed many interesting effects, different for different measures. The 
complete understanding of all possible sources of fluctuations and correla- 
tions in UrQMD would require much more detailed investigations and it is 
out of the scope of this publication. However, a rather consistent picture 
can be drawn from these basic tests: if one wants to study the exotic effects 
such as critical point, onset of deconfinement, etc. one should find good 
reference values of the studied measures. The UrQMD model suggests that 
the best reference values are those for negatively charged particles produced 
in very central A-\- A collisions, however even for these samples one should 
first try to better understand the origin of the rudimentary energy and/or 
system size dependence of A'''"'^ and T,^'^ observed in UrQMD. The simul- 
taneous measurement of all three quantities suggests that the system size 
(A) dependence of A'''"'^ and T,^'^ can be due to change of the single par- 
ticle spectrum mainly. The ^pj, values for negatively charged particles and 
for very central A + A are consistent with zero and do not depend on the 
mass of the colliding nuclei. For a sample of positively charged particles 
the reference values of all three measures may be masked by some trivial 
sources of fluctuations and correlations caused by the presence of protons in 
the sample. The fluctuation measures for a sample of all charged particles is 
always a non trivial combination of measures for positively and negatively 
charged particles. Moreover, correlations between positively and negatively 
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charged particles (i.e. due to particle decays) contribute here. 

Finally, one should mention that in most cases the restriction of the 
analysis to forward rapidity only reduces the values of fluctuation mea- 
sures. Therefore the NA61/SHINE experiment should do its best in order 
to measure event-by-event fluctuations in a rapidity window as wide as pos- 
sible. 
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